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We study the coupled-qubit oscillation driven by an oscillating field. When the period of the 
non-resonant mode is commensurate with that of the resonant mode of the Rabi oscillation, we 
show that the controUed-NOT (CNOT) gate operation can be demonstrated. For a weak coupling 
the CNOT gate operation is achievable by the commensurate oscillations, while for a sufficiently 
strong coupling it can be done for arbitrary parameter values. By finely tuning the amplitude of 
oscillating field it is shown that the high fidelity of the CNOT gate can be obtained for any fixed 
coupling strength and qubit energy gap in experiments. 
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I. INTRODUCTION 

The universal gate for quantum computing consists 
of the single qubit and the entangling two qubit oper- 
ations. Usually, electromagnetic oscillating fields such 
as microwave fields, laser pulses, and oscillating volt- 
ages have been used for single qubit operations. The 
atomic states of cavity-QED^ and ion-trap^ qubits are 
used as a natural basis of qubit states. In these cases 
the quantum Rabi-type oscillation can be analyzed in 
an approximation, called the rotating wave approxima- 
tion (RWA) . Similarly, the semi-classical Rabi-type oscil- 
lation of qubits of the artificial atomic states such as the 
superconducting charge qubit employed in the circuit- 
QED quantum computing'' and the flux qubit^ can also 
be analyzed in this approximation. 

We study the coherent two-qubit oscillation driven by 
an oscillating field. The two-qubit oscillation enables the 
two-qubit gate for the quantum computing. Among the 
two-qubit gates the CNOT gate is the most basic two- 
qubit operation for the universal gatei^ The CNOT gate 
operation was achieved in superconducting qubits with 
Ising-type interaction without driving oscillating field. ^ 
The oscillating-field-driven CNOT gate operation in su- 
perconducting qubits has recently been reported, but the 
fidelity is not so high due to the weak coupling strength 
between qubitsj^ In this study we propose a scheme for 
the CNOT gate operation between coupled qubits un- 
der an oscillating field for general Ising-type coupling 
strength. We discuss the CNOT gate operation for both 
the strong and weak coupling strength and show that the 
high fidelity of CNOT gate can be obtained even for a 
weak coupling. 

The CNOT gate uses the discriminating operations 
corresponding to different states of control qubit. De- 
pending on the control qubit states, the coupled-qubit 
state demonstrates a Rabi-type oscillation for the reso- 
nant oscillating field or non-resonant oscillation. During 
the TT rotation the target qubit state fiips to the other 
qubit state for a control qubit state, while for the other 
control qubit state the target qubit state goes back to the 
original state before reaching the other state. Though the 



latter oscillation is far from the Rabi-type oscillation, the 
commensurate oscillations give rise to the CNOT gate 
operation. By using the commensurate mode oscillations 
we obtain high fidelity for the CNOT gate operation. We 
show that for a weak coupling a high performance CNOT 
gate can be achieved by tuning the parameters, while for 
a sufficiently strong coupling the maximum fidelity can 
be obtained regardless of the values of system parame- 
ters. The maximum fidelity for a weak coupling can be 
obtained for any fixed coupling strength and qubit energy 
gap by finely tuning the amplitude of oscillating field in 
an experimental setup. This scheme of using the com- 
mensurate oscillation is quite general, so it is applicable 
to the natural atomic qubits as well as the solid state 
qubits. 



II. HAMILTONIAN OF COUPLED QUBITS 

The Hamiltonian of a two level system (qubit a) driven 
by an oscillating field with frequency uj is given by 



where 



El{K,t) 



g cos u>t, 



(1) 



(2) 



K is the external variable controlling the qubit energy 
levels, and a^^x are the Pauli matrices. Here g is the cou- 
pling strength between the qubit and the oscillating field 
which is proportional to the amplitude of the oscillat- 
ing field, iq is the tunnelling amplitude between different 
(pscudo-) spin states, and E°-{k) = E'^{k) — E'^{k). The 
qubit energy gap E°-{k) can be controlled by k, and for 
a particular value of kq the qubit can be brought to the 
degeneracy point, E'^[kq) = E'^{ko) = Eq. 

At this point the dominant energy scale is the tun- 
nelling energy and, if we introduce the coordinate 
transformation, |0) = (| |) -f | t))/V2 and |1) = (| | 
) — I t) )/ a/2, we have the Hamiltonian: 

H'^-£olO>(0|+£i11)(l|+.gcosc.i(|l)(OK|0)(l|), (3) 
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where S^^^^ 



Hence the resonant microwave 



with frequency to = 2t^ gives rise to the Rabi-type oscil- 
lation between the qubit states, |0) and |1). 

When two qubits (qubits a and b) are coupled, the 
Hamiltonian for coupled qubits in the basis of {| |), | f)} 
can be written as 



Here we introduce a coordinate transformation, 



H 



H"" (E)I + Ii 



tin 



with the Ising-type coupling strength J— 
1 



(4) 



(5) 



Here, Ess' are the energy levels of coupled qubits, 
EsAi^a, ^^b) = E'^^iKa) + EliKf,) ± J (+ for s = s', - for 
s = —s'), where Ka and Kb are the control variables for 
qubit a and 6, respectively. We here neglect the two-qubit 
tunneling term, t'^\s, s'){—s, —s'\, because it is negligibly 
small for usual parameter regimes. This term gives rise 
to the XY-type interaction which enables the SWAP gate 
operation rather than the CNOT gate. 
We can rewrite the Hamiltonian as 



H = J2Fss'i^ia,Hb, t)\s, S')(S, S'\ - i^|s, S'){-S, s'\ 

S.S' 



where — s is the opposite spin of s, and 



Euii^a, Hb) - 2gcosujt, 

Ei^iHa, Kb), 

E^^{Ka, Kb) + 2gcoaujt, 

E^liKa, Kb), 



(6) 



(7) 



Afterward, we will omit Ka, Kb in Ess'{Ka, Kb) for simplic- 
ity. 

To perform a CNOT gate operation, the system pa- 
rameters should be adjusted. We consider the situation 
that the external variables are adjusted in a way that 
(i) the energy gap /S.E = Ess' — E-ss' between different 
control qubit (qubit a) states is large and (ii) the tar- 
get qubit (qubit b) is at a degeneracy point, E\^\^ = E^-i;. 
This situation is usual for the CNOT gate operation of 
coupled qubits<^i^ In this case the tunneling process 
between different target qubit states takes a role, while 

between different control qubit states is negligible due 
to the large energy level difference /S.E. Consequently, 
the two-qubit Hamiltonian H becomes block-diagonal. 

Figure [Ija) shows the energy levels Ess' as a function 
of Kb, where we choose such that \Ess' ~ E^ss' \ ^ tq 
and thus i° can be negligible. In the figure there are two 
degeneracy points; lower degeneracy point where E^^ — 
E^-^ and upper degeneracy point where E^^ = E-^-^. By 
adjusting the variable Kb, the coupled-qubit states can 
be brought to one of these degeneracy points. Here the 
distance between these degeneracy points is related to 
the coupling strength between two qubitsi^i^ 



y = exp ( --o-y6'_^ je exp [-^(^v^t 



(8) 



with tan 61., = 2i^/|£;s; 



Es^\ in order to couple the os- 
cillating field with the off-diagonal elements of the qubit 
Hamiltonian. Then the Hamiltonian T-L = V^^HV be- 



comes 



n = [fpo(t)|pO)(pO| + £:pi(i)|pi)(pi| 

+apg COS ujt{\pO){pl\ + \pl){pO\)\, 



(9) 



where 



= C' - [(-!)" + {-^)''Pp\9COSiot (10) 

with uq — sin^j^, /Jq = cos0j_, ai — ainO^ and Pi = 
cos 9^. The two-qubit states \pp') are given as 

|pO) = cos(0,/2)|si)+sin(0,/2)|st>, 

\pl) = -sin(0,/2)|si)+cos(e,/2)|st>, (H) 



coupling 
strength 
11^-'"' 




FIG. 1: (Color online) (a) Energy levels Sppi of coupled qubits, 
where p, p' £ {0, 1}. E^s' with s, s' £ {^, t} are shown as thin 
dotted lines. The distance between two degeneracy points 
corresponds to the coupling strength between two qubits. (b) 
Occupation probabilities of \pp'} states during Rabi-type os- 
cillations at the lower degeneracy point where = E^^. 
Here we use the parameter values such that coupling strength 
J/h =0.6GHz, qubit energy gap a;o/27r=4GHz, and Rabi fre- 
quency Q.o/2ti = 600MHz. The initial state is chosen as 
^(0) = (|00) -I- |10))/\/2 and the CNOT gate is expected to 
be achieved at Clt = (odd) tt. 
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where s =1 (f) for p = 0(1). 

At the lower degeneracy point E^^_ = E^^ = Eq, we 
have the relations, 



(12) 



cO _ 
"--10(11) — 



E^^ + E^i I f E^^ — 



and 



9;=7r/2 and 61^ = tan"i(t^/2J) (13) 



Accordingly, the Schrodinger equation Hipit) = ih-^^{t) 
is written as ih§^(f>{t) = U(j){t) with % = U-^{t)nU{t) - 
ihU-^{t){dU{t)ldt) =no®'Hi, where 

o, _ ( 5icoswote"''^o*\ 

- ' cosc^ote^"«* ) ' ^^^^ 

uji — ujo gai coswote^*^^*-* 



gai cos cjote 



uji - Wo 



.(22) 



by using the relation of Eq. ([5]). Since the Hamiltonian 
7i is block-diagonal, we have 



-Ho = -^(|00)(00|-|01)(01|) 



Hi 



-gcoswt(|00)(01| + |01)(00|) 
Huji + 2(3igcoswt 



(|10)(10|-|11)(11|) 



-aigcosa;t(|10)(ll| + |11)(10|), 



(14) 
(15) 
(16) 



and ^(f) = Lo^t + 2(7/3i sinwoi/'^o- From this Hamiltonian 
we can obtain the two-qubit oscillation numerically. Also 
the dynamics can be analyzed in the RWA. 



III. ROTATING WAVE APPROXIMATION 

The RWA assumes near resonance w « cjq and weak 
coupling between a qubit and a oscillating field g/h <C 



10 



where 



ai 



= ^ine^ = ty^{2JY + [t\)\ (17) 



In the quantum Rabi oscillation for cavity-QED 
and ion-trap qubit, usually g/hujQ « 10^^ — 10^^. For 
the usual superconducting qubits the coupling strength 
g/hujQ ~ O(10~^) — — which is relatively strong, but we 
find that the RWA gives accurate results consistent with 
our numerical calculation. 

The off-diagonal element of 'Ho(i) is written as 



/3i = cos0t = 2J/^(2J)2 + (t^)2 

Here, aig corresponds to the transition frequency be- 
tween 1 10) and 1 11) states, while /3i.g term induces un- 
necessary complicate oscillations. In Fig. [T]the energy 
gaps are given as 



-Hi 



= 1(1 
12 2^ 

^gai 
12 2 



^ V ^0 



(23) 
,(24) 



where Jn{x) is the Bessel function of the first kind. In 
the usual RWA, the fast oscillating term e^'""* in [Ho] 12 
is neglected. Similarly we here neglect e™'^"* (n 0) in 
1 = - £?o = 2^[{E,^~Eu)/2f + {tb^y, [^i]i2, resulting in 



n^o = fo°i-C = 2t^ (18) 



where Hujq is the qubit energy gap, and hwi depends on 
the qubit coupling strength through the relation 



-TV RWA 



UJ-. = UJr, + \ — 



(19) 



We here consider the case that the oscillating field is 
resonant with the energy gap wq between the states |00) 
and |01), i.e. oj = ujq, at the degeneracy point En — 
E^^. Then this Hamiltonian describes the usual Rabi- 
type oscillation between the states |00) and |01) with 
the Rabi frequency « ff/^, while the evolution of 
the states |10) and |11) is far from the Rabi oscillation, 
since the energy level difference is not resonant with the 
oscillating field frequency, uji =/= lu, for a finite coupling 
strength J. 

We again introduce a rotating coordinate such as 
tpit) = U{t)(j){t), where 



RWA_ 



.9/2 
5/2 
h{uji ~ Wo) g'/2 

g' /2 h{uji - ujq) 



(25) 
(26) 





i 






exp 




® exp 





2g/?i 

UJQ 



sincjo^ c 



(20) 



where g' = gai [Jq {2g(3i/fu^o) + J-2 i2g/3i/hwo)] . 

Hence the Hamiffonian H^^A ^ (l/2)ga.^ = {9/^)8^ 
describes the Rabi oscillation with the Rabi frequency 
JIr = fio = g/h, while the Hamiltonian lif^^ shows 
a non-resonant oscillation with the oscillating frequency 

ill = \J (t^i — ^of' + {g' I^Y- From the relation of Eq. 
([T9| we see that the behavior of this non-resonant oscil- 
lation depends on the coupling strength J as well as wq 
and g. 

In Fig. [Ijb) we show the resonant and non-resonant 
oscillations, when Kf, is adjusted to the lower degener- 
acy point where E^ ~ E^. Then a microwave with 
resonant frequency w = wq gives rise to the Rabi oscil- 
lation between two states |00) and |01), while the states 
1 10) and 1 11) experience a non- resonant oscillation. The 
controlled-NOT gate operation requires that the target 
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qubit flips for a specific state of control qubit such that 
|00) ^ |01) while |10) ^ |10). However, for example, at 
ilt — TT in Fig. [lib) the states |11) and |10) also evolves 
during the transition from |00) to |01). Thus we cannot 
expect a good CNOT gate operation in this case. 

Although, for coi different from the resonant value of 
Lo = luq, the oscillation is not a Rabi oscillation, the os- 
cillation period can be an even integer multiple of that of 
the resonant Rabi oscillation mode for a specific values 
for parameters g, J, and ujq which correspond to the os- 
cillating field amplitude, the coupling strength between 
qubits, and the qubit energy gap, respectively. The con- 
dition for this commensurate periods is given by 



„ 27r 27r 
2n— = — , 



(27) 



as we can see in Fig. [5] This condition determines the 
value of g for given values of (wq, J)- 

From Eq. (|27p the value of g for fidelity maxima can 
be expressed as 



h^LJl - Wo) 



(28) 



The argument of the Bessel function is written as 
2gl3i/hujo = (2g/fta;o)4J/v/(M0^TT6J2. For x ^ 0, 
the Bessel functions Jo (a;) and J-2{x) approach 1 and 0, 
respectively. Thus, for small J and large ojq the expres- 
sion of g in Eq. (j28|) can be approximated as 



1 



9 ' 



hujQ 



(29) 



using Eq. ([HI). These expressions of Eqs. ([28|) and ([29]) 
provide the value for g for the fidelity maxima with given 
values of J and luq. 
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FIG. 3: (Color online) (a) Fidelity F for CNOT gate as a func- 
tion of g at Q,ot — n. The main (dark arrow) and subsidiary 
(light arrow) maxima correspond to the commensurate modes 
in Fig. [2fa) and (b), respectively. Here we set J/h — 0.5GHz 
and UJQ = 4GHz. Inset shows the fidelity envelops as a 
function of cjq for g/h — 0.265GIIz and J/h — O.SGHz. 
(b) Fidelity F as a function of J for g/h = 0.265GIIz and 
Ldo ~ 4GIIz. Inset shows F for different g's. 




FIG. 2: (Color online) (a) Commensurate oscillations of oc- 
cupation probability of coupled-qubit states with the initial 
state, |t/j(0)) = (|00) + |10))/^/2 for g/h = 0.265GHz. The 
non-resonant oscillation modes (Pio and Pn) are commensu- 
rate with the resonant modes (Poo and Poi). At Qt = (odd)7r, 
Pio and Pii recover their initial values, thus the CNOT gate 
operation is achieved. Here Qo = g/h, J/h = 0.5GHz, and 
loq/2-k =4.0GHz. (b) Higher order commensurate modes for 
smaller g/h — 0.122GHz with the same J and luq- 



IV. CNOT GATE OPERATION USING 
COMMENSURATE MODES 

The scheme for CNOT gate operation in this study 
uses the non-Rabi oscillations for |10) and |11) states 
which are commensurate with the Rabi oscillation for |00) 
and |01) states. In Fig. [2|we display the numerical results 
obtained from the Hamiltonian in Eqs. (|21l) and (|22p . 
which show such commensurate mode oscillations. The 
initial state, |V'(0)) = (|00) -I- |10))/V2, is driven by an 
oscillating field with the resonant frequency oj — luq < uji. 

In experimental situations usually the coupling J and 
the qubit energy gap wq are set to be fixed. Thus the con- 
trol of oscillating field amplitude g with fixed J and ojq 
will be more desirable. For any given pair of {loq , J) one 
can find a commensurate oscillation by finely tuning the 
value of g according to Eq. (1^51) . By varying g with fixed 
ujQ and J, we were able to find a commensurate oscilla- 
tion mode numerically; the oscillation period of the |00) 
and |01) states is twice of that of the |10) and |11) states 
[Fig. [H^a)], which corresponds to n = 1 in Eq. ([?7| . As 
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FIG. 4: (Color online) Rabi-type oscillations of occupation 
probabilities of \pp') states for strongly coupled qubits with 
the initial state '0(0) = (|00) + |10))/\/2- Here the parameters 
are J jh = 5GHz, cjo/27r=4GHz, and Q.qI'I-k = 600MHz at the 
degeneracy point where E^^ — E^^ in Fig. [T][a). 



g decreases further, another commensurate mode with a 
shorter period appears in Fig. ^h) {n = 2). Actually 
we have found a series of commensurate modes as g de- 
creases. The values of g obtained numerically coincide 
well with those from the RWA in Eq. (|28p as shown in 
Table HI 

The CNOT gate operation is done when the occupation 
probability Pqo (^oi) is reversed perfectly from 0.5 (0) to 
(0.5) at fit — (odd) it. At the same time, we can 
observe that the probabilities Pio and Pu recover their 
initial values 0.5 and 0, respectively. As a result, the 
CNOT operation is realized by using these commensurate 
oscillations. 

Let us consider a concrete example for comprehen- 
sive understanding. For superconducting flux qubits;^!^ 
g = mB is the coupling between the amplitude B of the 
magnetic microwave field and the magnetic moment m, 
induced by the circulating current, of the qubit loop. In 
order to adjust the value of g, actually we need to vary the 
microwave amplitude B, because the qubit magnetic mo- 
ment is fixed at a specified degeneracy point. The Rabi- 
type oscillation occurs between the transformed states 
|0) = (U) + I t))/V2 and |1) = (U) - I t))/V2. The 
states of qubits can be detected by shifting the magnetic 
pulse adiabaticallyi^. Since these qubit states are the su- 
perposition of the clockwise and counterclockwise current 
states, I I) and | f), the averaged current of qubit states 
vanishes at the degeneracy point in Fig. [T{a). Thus, one 
can apply a finite dc magnetic pulse to shift the qubits 
slightly away from the degeneracy point to detect the 
qubit current states. 

We now discuss the performance of CNOT gate oper- 
ation. The fidelity for CNOT gate operation is given by 
F{t) = Tr(M(t)McNOT)/4*^ where Mcnot is the matrix 
for the perfect CNOT operation and M{t) is the truth 



table amplitude at time t. In Fig. [3] we plot the fidelity 
F by varying g, luq or J at ftot = tt. In Fig. the 
main and subsidiary maxima correspond to the CNOT 
operation in Fig. [IJa) and (b), respectively. As shown in 
Eq. (I28p the fidelity maxima are determined by the three 
parameters, g, ujq or J. The series of maxima correspond 
to the different n in Eq. (^5)) . In the inset of Fig. (HJa) 
we also show the fidelity as a function of ujq. 

An interesting behavior of fidelity maxima is shown 
in Fig. ^h) as a function of J, where the fidelity ap- 
proaches the maximum as the coupling strength J in- 
creases. In the inset we show the oscillations for various 
parameters converge to 1, which implies that for suffi- 
ciently strong coupling maximum fidelity for the CNOT 
gate is achievable regardless of the values of other param- 
eters. This is because, for a sufficiently strong coupling 
J 3> wq, thus ai — sin^i « 0, the off-diagonal terms in 
Eq. (|16p which induce the oscillation between two states 
with p — 1 are vanishing and thus the occupation proba- 
bilities of the 1 10) and |11) states are not changed. As a 
result, the states |10) and |11) preserve their initial occu- 
pation probabilities, while the states |00) and |01) expe- 
rience a Rabi-type oscillation. In this limit, the CNOT 
gate operation can be achieved with arbitrary parameter 
values. 

In Fig. Slwe show the Rabi-type oscillation for strongly 
coupled qubits. While the Pqo (Pqi) is reversed from 0.5 
(0) to (0.5) at iU — (odd) tt, we can observe that the 
probabilities Piq and Pu remain their initial values 0.5 
and 0, respectively. In this case the parameters need not 
satisfy the commensurate condition of Eq. (P7)) for the 
CNOT gate operation. 

In fact, however, it is not so easy to obtain a sufficiently 
strong coupling between qubits in experiment. Instead, 
we can achieve a high fidelity of CNOT gate by choosing 
parameters satisfying Eq. (j28p . In real experiments we 
can control the amplitude of the oscillating field while 
the coupling J and energy gap luq are fixed. In Table 
H] we compare the numerical results for g at the main 




FIG. 5: (Color onUne) The errors of fidelity, SF ^ I ~ F, at 
the main peak for various values of J and ujo at f2oi = tt. 
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0.011 


0.100 
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RWA 


0.011 


0.100 


0.264 


0.484 
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g/h 




numerical 


0.023 


0.185 


0.448 


X 


X 




=2) 


RWA 


0.023 


0.184 


0.443 


0.752 


1.090 



TABLE I: The values of g/h for the main fidehty maxima 
(n = 1) obtained from numerical calculation and from the 
RWA of Eq. (|28p for various coupling J and qubit energy gap 
ojQ. For small uio and large J the oscillations are far from the 
Rabi oscillation. Here, the unit of all numbers is GHz. 

fidelity maximum points {n—1) with those obtained from 
Eq. ([28l) in the RWA, which shows that two values fit 
well with each other. The approximate values are a little 
smaller than the numerical values and the small deviation 
tends to increase as J increases and wq decreases. This 
can be understood from Eq. (|29|) where g increases as J 
increases and ujq decreases. Hence for small J and large 
Wo the RWA and the numerical calculation coincide with 
each other, because the RWA works well in the regime 
g/h-^ wq. For the parameters far away from this regime 
the two-qubit oscillation deviates seriously from the Rabi 
oscillation and thus the CNOT gate operation cannot be 
performed, except the strong coupling limit discussed. 

Though the value of fidelity at the main peak in Fig. [3] 
is close to the maximum value of 1, it has small deviation. 



5F = 1 — F . In Fig. [5] we plot the fidelity error 5F 
for various values of [ujq^J) at 0.^1 = tt. For large loq 
and small J the fidelity error is vanishing, and the Rabi 
oscillation by the Hamiltonian T-L in Eqs. (I?T|) and ([2^ 
is close to 5F sa 10^^ for the fault-tolerant quantum 
computing. Hence for a weak coupling as well as a strong 
coupling we can achieve high performance CNOT gate 
operation. 



V. SUMMARY 

The commensurate oscillations of resonant and non- 
resonant modes enable the high fidelity CNOT gate op- 
eration by finely tuning the oscillating field amplitude 
for any given values of qubit energy gap and coupling 
strength between qubits. While for a sufficiently strong 
coupling the CNOT gate can be achieved for any given 
parameter values, for a weak coupling a relation between 
the parameters should be satisfied for the fidelity max- 
ima. For a sufficiently weak coupling compared to the 
qubit energy gap, J/huJo ^ 1, we have ai ~ 1 and 
/?! ~ 0, resulting in the expression for g in Eq. (f^ . 
For J/hbjQ ^ 1, Eq. ((29|) immediately gives rise to the 
relation g/ J <ti 1 and thus g/hujQ <C 1 after some manip- 
ulation. This means that for a weak coupling J/Hujq <C 1 
the numerical results are well fit with the RWA as shown 
in Table m because the RWA is good for g/hojQ <C 1. As 
a result, the high performance CNOT gate operation can 
be achieved as shown in Fig. [5l 
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